RBRC-856 

Photon decay in strong magnetic field in heavy-ion collisions 

Kirill Tuchin a ' fe 

a Department of Physics and Astronomy, 
Iowa State University, Ames, IA 50011 
b RIKEN BNL Research Center, 
Upton, NY 11973-5000 



(Dated: February 18, 2011) 

We calculate the photon pair production rate in strong magnetic field created in off-central 
heavy-ion collisions. Photon decay leads to depletion of the photon yield by a few percent 
at RHIC and by as much as 20% at the LHC. It also generates a substantial azimuthal 
asymmetry ( "elliptic flow" ) of the final photon distribution. We estimate v 2 ~ 2% at RHIC 
and v 2 ~ 14% at LHC. Photon decay measurements is an important tool for studying the 
magnetic fields in early stages of heavy-ion collisions. 

Ultra-relativistic heavy ions colliding at finite impact parameter possibly create a super-critical 
magnetic field B. According to the estimates in [1, 2], the strength of this field at yfs = 200 GeV 
is about eB ~ rri^/h, while the critical field is eB c = m\j%. Thus, magnetic field created in 
heavy-ion collisions is by many orders of magnitude stronger than any field that has been created 
using the state-of-the-art lasers (see e.g. [3]). Possible existence of such fields opens a new avenue 
for studying the high intensity regime of QED. 

Various QED processes in external magnetic field strongly depend on the time-dependence of 
that field. Recently we argued [4] that the magnetic field is approximately stationary during 
the life-time of the quark-gluon plasma (QGP) that is formed shortly after the collision. Indeed, 
phenomenological models describing evolution of QGP indicate that the thermalized medium is 
formed almost immediately after the collision (after ~ 0.5 fm [5, 6]) when the magnetic field is 
near the maximum of its strength. As the heavy-ion remnants recede from the collision point, the 
magnetic field tends to rapidly decrease with time. This induces circular currents in the QGP that, 
by the Faraday law, produce an induced magnetic field in the direction of the external field. Thus, 
the relaxation process of the external field slows down. The characteristic relaxation time is [4] 



where R is the QGP size and a is its electric conductivity. Lattice calculations show that the 
electric conductivity is high even at temperatures close to T c [7]. In [4] we used the lattice data 
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of Ref. [7] to estimate the relaxation time as r « 2.2 fm (T/T c ) 2 . This number is even larger if 
the effect of magnetic field on the electric conductivity is taken into account [8]. It implies that 
external magnetic field is a slowly varying function of time during the entire QGP life-time. Of 
course, once the plasma cools down to the critical temperature and undergoes the phase transition 
to the hadronic gas, the conductivity becomes very small and the magnetic field cannot be sustained 
anymore. 

In [4] we discussed the properties of the synchrotron radiation of gluons by fast quarks and 
argued that it has a significant phenomenological implications. Indeed, the corresponding energy 
loss in magnetic field is comparable to that sustained by the fast quark in hot nuclear medium. 
The azimuthally asymmetric form of the energy loss contributes to the 'elliptic flow' phenomenon 
observed at RHIC. In this letter we consider a cross-channel process - pair-production by photon 
in external magnetic field. Specifically, we are interested to determine photon decay rate w in the 
process "yB — > ffB, where / stands for a charged fermion, as a function of photon's transverse 
momentum kx, rapidity r\ and azimuthal angle ip. Origin of these photons in heavy- ion collisions 
will not be of interest in this paper. 

Characteristic frequency of a fermion of species a of mass m a and charge z a e (e is the absolute 
value of electron charge) moving in external magnetic field B (in a plane perpendicular to the field 
direction) is 



where e is the fermion energy. Here - in the spirit of the adiabatic approximation - B is a slow 
function of time. Calculation of the photon decay probability significantly simplifies if motion of 
electron is quasi-classical, i.e. quantization of fermion motion in the magnetic field can be neglected. 
This condition is fulfilled if Hub <C e. This implies that 



For RHIC it is equivalent to e > m n , for LHC e S> 4m n . 

Photon decay rate was calculated in [9] and, using the quasi-classical method, in [10]. It reads 



e > VzeB . 



(3) 



w = — 
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where summation is over fermion species and the invariant parameter k is defined as 





( V) 2 = 




(k x B) 2 , 
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with the initial photon 4-momentum k^ = (Huj, k). In heavy-ion collisions the vector of magnetic 
field B is orthogonal to the "reaction plane" , which is spanned by the impact parameter b and the 
collision axis z. We define the polar angle 9 with respect to the z-axis and azimuthal angle ip with 
respect to the reaction plane. In this notation, B = By and k = k z z + k± (x cos ip + y sin ip), where 
k± = \k\sin6 = friosinO. Thus, (B x k) 2 = B 2 (k 2 + k^cos 2 tp). Introducing rapidity rj as usual 
hu = k± cosh rj and k z = k± sinh rj we can write 

h(z a eB) 



rrvi 



k± y sinh rj + cos 2 ip . 



(6) 



In Fig. 1 we plotted the photon decay rate (4) for RHIC and LHC. The survival probability of 
photons in magnetic field is P = 1 — wAt, where At is the time spent by a photon in plasma. We 
can see that for At = 10 fm the photon survives with probability Prhic ~ 97% at RHIC, while 
at LHC Plhc ~ 80%. Such strong depletion can certainly be observed in heavy-ion collisions at 
LHC. 
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FIG. 1: Decay rate of photons moving in reaction plane in magnetic field as a function of transverse 
momentum kr'- (a) at RHIC, (b) at LHC. Broken lines from bottom to top give contributions of 7 — ¥ dd 7 
7 — > uu, 7 — > ji + ji~ and 7 — > e + e~ channels. Upper solid line is the total rate. 



Azimuthal distribution of the decay rate of photons at LHC is azimuthally asymmetric as can be 
seen in Fig. 2. The strongest suppression is in the B field direction, i.e. in the direction orthogonal 
to the reaction plane. At rj > 1 the ip dependence of K a is very weak which is reflected in nearly 
symmetric azimuthal shape of the dashed line in Fig. 2. 

To quantify the azimuthal asymmetry it is customary to expand the decay rate in Fourier series 
with respect to the azimuthal angle. Noting that w is an even function of tp we have 



w 



j 00 f'7T 

(f) = 2 W ° + ^2 w n cos(nip) , w n = - 

n=l 71 ^~ 7r 



w{ip) cos(mp) dip . 
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FIG. 2: Azimuthal distribution of the decay rate of photons at different rapidities at LHC. Only contribution 
of the 7 — > e + e~ channel is shown. 



In strong fields x a S> 1. For example, for 7 — > fi + fi at RHIC at <p = rj = and kx = 1 GeV we 
get x M = 19. Therefore, we can expand the rate (4) at large x a as [9] 

A 



w 



3 1/6 5r 2 (|) v a cm effz Q 3 _ ^ 

2 4/3 77r l/2 r (7) ^^1/3 (sinh^ + COsV) 1 / 6 ' 



(8) 



At rj = the Fourier coefficients w n can be calculated analytically using formula 3.631.9 of [11] 

3 2 1 /3,4 
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, ^2fc+i = 0, fc = 0, 1,2, ... , 



(9) 



where B is the Euler's Beta- function and A is defined in (8). Substituting these expressions into 
(7) we find 
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cos(2k(p) 



The first few terms in this expansion read 



1 / 2 . 14 

w = -w I 1 - - cos(20) + — cos(40) 



(10) 
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What is measured experimentally is not the decay rate, but rather the photon spectrum. This 
spectrum is modified by the survival probability P which is obviously azimuthally asymmetric. To 
quantify this asymmetry we write using (7) 

l-P2w 2k 

V 2 k COS(Z(fK) I , V 2 k = - 

1,-1 J 
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where P = (1 — wAt)^ = 1 — w^At is the survival probability averaged over the azimuthal angle. 
Since woAt <C 1, as can be seen in Fig. 1, we can estimate using (8) and (9) 

^ 2w 2k 2 W2fc/W 5 6 2 / 3 r(§) ^ a cm (eB) 2 / 3 z S J 3 
v 2 k~ w At = At —y T u/3 • (13 

Wo Wo ITT (Kt) I 

In particular, the "elliptic flow" coefficient is 

For example, at kx = 1 GeV and At ~ 10 fm/c one expects v 2 — 2% at RHIC and v 2 — 14% 
at LHC only due to the presence of the magnetic field. We see that decay of photons in external 
magnetic field significantly contributes to the photon asymmetry in heavy-ion collisions along with 
other possible effects [12-17]. 

In summary, we calculated photon pair-production rate in external magnetic field created in 
off-central heavy-ion collisions. Photon decay leads to depletion of the photon yield by a few 
percent at RHIC and by as much as 20% at the LHC. The decay rate depends on the rapidity 
and azimuthal angle. At mid-rapidity the azimuthal asymmetry of the decay rate translates into 
asymmetric photon yield and contributes to the "elliptic flow". Let us also note that photons 
polarized parallel to the field are 3/2 times more likely to decay than those polarized transversely 
[9]. Therefore, polarization of the final photon spectrum perpendicular to the field is a signature 
of existence of the strong magnetic field. Finally, photon decay leads to enhancement of dilepton 
yield, which will be addressed in a separate publication. 
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